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We review the main features of the generalized Chaplygin gas (GCG) proposal for 
unification of dark energy and dark matter and discuss how it admits an unique 
decomposition into dark energy and dark matter components once phantom-like 
, dark energy is excluded. In the context of this approach we consider structure 

Q ■ formation and show that unphysical oscillations or blow-up in the matter power 

I ' spectrum are not present. Moreover, we demonstrate that the dominance of dark 

O ' energy occurs about the time when energy density fluctuations start evolving away 

from the linear regime. 



1. Introduction 



The GCG model 1 ' 2 is an interesting alternative to more conventional ap- 
proaches for explaining the observed accelerated expansion of the Universe 
such as a cosmological constant 3 or quintessence 4 . It is worth remarking 
that quintessence is related to the idea that the cosmological term could 
evolve 5 and with attempts to tackle the cosmological constant problem. 

In the GCG approach one considers an exotic equation of state to de- 
scribe the background fluid: 

Pch = 4" ' W 

Pch 

where A and a are positive constants. The case a = 1 corresponds to the 
Chaplygin gas. In most phenomenological studies the range < a < 1 



"Talk presented at the V New Worlds in Astroparticle Physics Conference, Faro, Portu- 
gal, 8-10 January 2005. 
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is considered. Within the framework of Friedmann-Robertson- Walker cos- 
mology, this equation of state leads, after being inserted into the relativistic 
energy conservation equation, to an evolution of the energy density as 2 

i 

l+a 

(2) 

where a is the scale-factor of the Universe and B a positive integration 
constant. From this result, one can understand a striking property of the 
GCG: at early times the energy density behaves as matter while at late 
times it behaves like a cosmological constant. This dual role is what es- 
sentially allows for the interpretation of the GCG model as an entangled 
mixture of dark matter and dark energy. 

The GCG model has been successfully confronted with different classes 
of phenomenological tests: high precision Cosmic Microwave Background 
Radiation data 6 , supernova data 7 , and gravitational lensing 8 . More re- 
cently, it has been shown using the latest supernova data 9 , that the GCG 
model is degenerate with a dark energy model with a phantom-like equa- 
tion of state 10 ' 11 . Furthermore, it can be shown that this does not require 
invoking the unphysical condition of violating the dominant energy condi- 
tion and does not lead to the big rip singularity in future 10 . It is a feature 
of GCG model, that it can mimic a phantom-like equation of state, but 
without any kind of pathologies as asymptotically the GCG approaches to 
a well-behaved de-Sitter universe. Structure formation has been studied in 
Refs. [2,12]. In Ref. [13], the results of the various phenomenological tests 
on the GCG model are summarized. 

Despite these pleasing performance concerns about such an unified 
model were raised in the context of structure formation. Indeed, it has 
been pointed out that one should expect unphysical oscillations or even 
an exponential blow-up in the matter power spectrum at present 14 . This 
difficult arises from the behaviour of the sound velocity through the GCG. 
Although, at early times, the GCG behaves like dark matter and its sound 
velocity is vanishingly small as one approaches the present, the GCG starts 
behaving like dark energy with a substantial negative pressure yielding a 
large sound velocity which, in turn, produces oscillations or blow-up in the 
power spectrum. In any unified approach this is inevitable unless the dark 
matter and the dark energy components of the fluid can be properly iden- 
tified. These components are, of course, interacting as both are entangled 
within a single fluid. However, it can be shown that the GCG is a unique 



Pch = 



A + 



B 



a 3(l+c) 



February 2, 2008 17:35 Proceedings Trim Size: 9in x 6in 



faro05 



3 

mixture of interacting dark matter and a cosmological constant-like dark 
energy, once one excludes the possibility of phantom- type dark energy 15 . 
It can be shown that due to the interaction between the components, there 
is a flow of energy from dark matter to dark energy. This energy transfer is 
vanishingly small until recent past, resulting in a negligible contribution at 
the time of gravitational collapse (z c ~ 10). This feature makes the model 
indistinguishable from a CDM dominated Universe till recent past. Subse- 
quently, just before present [z ~ 2), the interaction starts to grow yielding 
a large energy transfer from dark matter to dark energy, which leads to the 
dominance of the latter at present. Moreover, it is shown that the epoch of 
dark energy dominance occurs when dark matter perturbations start devi- 
ating from its linear behaviour and that the Newtonian equations for small 
scale perturbations for dark matter do not involve any fc-dependent term. 
Thus, neither oscillations nor blow-up in the power spectrum do develop. 



2. Decomposition of the GCG fluid 

In Ref. [2], it is shown that the GCG can be described through a complex 
scalar field whose Lagrangian density can be written as a generalized Born- 
Infeld theory: 



L G bi = 



i - {g^e^e, v )- 



(3) 



which reduces into the Born-Infeld Lagrangian density for a = 1. The field 
9 corresponds to the phase of the complex scalar field 2 . 

Let us now consider the decomposition of the GCG into components. 
Introducing the redshift dependence and using Eqs. (1) and (2), the pres- 
sure is given by 



Pch = -- 



[A + B(l + z) 3 ( 1 + Q )] !+« 
while the total energy density can be written as 



(4) 



Pch = \A + B(l + zf^] 1+ 



(5) 



where the present value of the scale- factor, oo, has been set to 1. 

We decompose the energy density into a pressure-less dark matter com- 
ponent, Pdmi and a dark energy component, px, with an equation of state 
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Px — wxpx] hence the equation of state parameter of the GCG can be 
written as 



Pch Px w x px 



Pch Pdm + PX Pdm + PX 

Therefore, from Eqs. (4), (5) and (6), one obtains for px 



(6) 



Pdm /_n 

PX ~ 1 + W X [l + f(l + s)3(l+a)] ' ^ > 

From this equation one can see that requiring that px > leads to the 
constraint wx < for early times (z 3> 1) and wx < — 1 for the future 
(z = —1). Thus, one can conclude that Wx < —1 for the entire history of 
the Universe. The case wx < — 1 corresponds to the so-called phantom-like 
dark energy, which violates the dominant-energy condition and leads to an 
ill defined sound velocity (see however Ref. [10]). Excluding this possibility, 
then the energy density can be uniquely split as 



where 



P = Pdm + PA (8) 



[A + B(l + z) 3 ( 1 +«)] i+o 



and 



^4 

Pa = -Pa = = «r , (10) 

[A + B{1 + z)3( 1 +«)] !+« 

from which one finds the scaling behaviour of the energy densities 

&H = 4(l + z) 3 ( 1+Q ) . (11) 
PA A 

In what follows we express parameters A and B in terms of cosmological 
observables. From Eqs. (9) and (10), it implies that 



l 

PchO = PdmO + Pao = {A + B) i+« , (12) 
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Figure 1. Density parameters Qdm an d £Ia and SI;, as a function of redshift. It is 
assumed that Q b0 = 0.05, f^dmO = 0.25, Hao = 0-7 and a = 0.2. 



where p c ho, Pdma and pao are the present values of p c u, p m and p\, respec- 
tively. Constants A and B can then be written as a function of p c ho 



A = pao p\ 



chO i 



B — PdmO PchO 



(13) 



It is also interesting to express A and B in terms of f^mo, £^ao, the 
present values of the fractional energy densities Qdm(A) — Pm(A)/Pc where 
p c is the critical energy density, p c — 3H 2 /8irG. From the Fricdmann 
equation 



3H 2 = 8ttG 



i_ 

A + B(l + z) 3{1+a) 1+a +8irG Pb0 (l + zf 



(14) 



where pto is the baryon energy density at present, one obtains 



A ~ n A0 P ( c l +a) , B ~ n dm0 p { ^ +a) . (15) 
Therefore, with the present value of the Hubble parameter, Hq: 
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H 2 = Hi 



i_ 

^ao + n dm0 (l + ^) 3(1+a) 1+a + "feo(l + z) 



one can write the fractional energy densities Q d m , ^A and as 



(16) 



[n AO + ^ dm o(i + z)W] a/(1+a) x 

^A0 



» fc0 (l + z) 3 
X 



where 



n A0 + a„ (i + ^) 3(1+a) 



1/(1+") 



+ Q b0 (l + z) s 



(17) 



(18) 



(19) 



(20) 



Finally, as fidmo and il\o are order one quantities, one can easily see 
that at the time of nucleosynthesis, J7a is negligibly small, and hence the 
model is not in conflict with known processes at nucleosynthesis. 

It is important to realize that there is an explicit interaction between 
dark matter and dark energy. This can be understood from the energy 
conservation equation, which in terms of the components can be written as 



Pdm + 3Hp dm = -p A . (21) 

Thus, the evolution of dark energy and dark matter are coupled so that 
energy is exchanged between these components (see Refs. [16,17] for earlier 
work on the interaction between dark matter and dark energy). One can 
see from Figure 1, that until z ~ 2, there is essentially no exchange of 
energy and the A term is vanishingly small. However, around z ~ 2, the 
interaction starts to increase, resulting in a substantial growth of the dark 
energy term at the expense of the dark matter energy. Thus, by around 
z ~ 0.2, dark energy starts dominating the energy content of Universe. Of 
course, these redshift values are a dependent and, in Figure 1, a = 0.2 has 
been chosen. Nevertheless, the main conclusion is that in this unified model, 
the interaction between dark matter and dark energy is vanishing small for 
almost the entire history of the Universe making it indistinguishable from 
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the CDM model. As can be clearly seen, the energy transfer has started in 
the recent past resulting in a significant energy transfer from dark matter to 
the A-like dark energy. In the next section we show that this energy transfer 
epoch is the one when dark matter perturbations start departing from its 
linear behaviour. But before that notice also that Eq. (21) expresses the 
energy conservation for the background fluid, which is reminiscent of earlier 
work on varying A cosmology 5, 18,19 where the cosmological term decays into 
matter particles. In here, we have the opposite, as a is always positive, 
hence the energy transfer is from dark matter to dark energy. This is 
responsible for the late time dominance of the latter and ultimately to the 
observed accelerated expansion of the Universe. 



3. Structure Formation 

Aiming to study structure formation, it is interesting to write the 0-0 com- 
ponent of Einstein's equation as 

'iH 2 = 8irG(pdm + Pb) + A , (22) 

where A is given by 



A = 8ttGpa • (23) 

We address now the issue of energy density perturbations. We first write 
the Newtonian equations for a pressure-less fluid with background density 
Pdm and density contrast 5d m , with a source term due to the energy transfer 
from dark matter to dark energy. Assuming that both, the density contrast 
5dm and the peculiar velocity v are small, that is 5dm << 1 and v « u, 
where u is the velocity of a fluid element, one can write the Eulcr, the 
continuity and the Poisson's equations in the co-moving frame 19 : 



dv a V$ 

ax + — + -v = 

at a a 



V • v 



r V 2 $ 



85, 



dm 



(//;/ 



dt 

^TlGpdm{l 



Pdn 



Sdr, 



A 



(24) 
(25) 
(26) 



where $ is the gravitational potential, and W is the source term in the 
continuity equation due to the energy transfer between dark matter and 



February 2 



faro05 



10 - 



0.1 ^ 



0.01 



0.001 




0.001 



Figure 2. Density profile <5,j m as function of scale factor. The solid, dotted, dashed 
and dash-dot lines correspond to a = 0, 0.2, 0.4, 0.6, respectively. It is assumed that 
n b0 = 0.05, n dm0 = 0.25 and n A0 = 0.7. 



the cosmological constant-type dark energy. The co-moving coordinate x 
is related to the proper coordinate rbyr = ax. In here, 



* = — 



8ttG 



:A 



(27) 



One expects a perturbation also in the A term. However, it can be seen 
from the Euler equation, for a fluid with an equation state of the form 
P = wp, 



(w + l)p + v ■ Vv^ + wVp + (w + l)pV$ = (28) 

so that, for w = — 1, it follows that Vp = 0, from which implies that this 
cosmological constant like component is always homogeneous. We should 
mention that the Euler Eqs. (26) and (28) can have an extra term in the 
r.h.s. if the velocity of the created A-like particle has a different velocity 
from the decaying dark matter particle 19 . In this case, the A-like dark 
energy can have spatial variations which can be neglected for the Newtonian 
treatment. However, in our case, we are considering only the situation 
where both the decaying and created particles have the same velocity. 
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Figure 3. The growth factor m(y) as a function of scale factor a. The solid, dotted, 
dashed and dash-dot lines correspond to a = 0,0.2,0.4,0.6, respectively. It is assumed 
that Q b0 = 0.05, n dm0 = 0.25 and Q A0 = 0.7. 



From the divergence of Eq. (24) and using Eqs. (25) and (26), one 
obtains the small scale linear perturbation equation for the dark matter in 
the Newtonian limit: 



d 2 8 dri 
dt 2 



+ 



~a * 
2- + 


dS dm 


a pdm_ 


dt 



AnGpdm - 2- 



a Pdm 



d_ 

di 



Pdm 



Jdm 



(29) 

One sees that, if \I/ = 0, that is in the absence of energy transfer, one 
recovers the standard equation for the dark matter perturbation in the 
ACDM case. One can verify that this occurs for a = 0. It can also be seen 
from the above equation that there is no scale dependent term to drive 
oscillations or to cause any blow up in the power spectrum. 

We turn now to the evolution for the baryon perturbations in the New- 
tonian limit when the scales are inside the horizon. Given that our purpose 
is to consider the period after decoupling, the baryons are no longer coupled 
to photons and one can effectively consider baryons as a pressure-less fluid 
like the dark matter as there is no significant pressure due to Thompson 
scattering. We assume that there is no interaction between dark energy 
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and baryons, which means that the Equivalence Principle is violated as, on 
its turn, dark matter and dark matter are strongly coupled. Given that 
it is parameter a that controls this interaction (a = means there is no 
interaction), it is a measure of the violation of the Equivalence Principle. 
One can also see from the behaviour of that this violation also starts 
rather late in the history of the Universe. In the Newtonian limit, the evo- 
lution of the baryon perturbation after decoupling for scales well inside the 
horizon is similar to the one for dark matter however, as described earlier, 
the source term is absent as there is no energy transfer to or from baryons. 
Thus, the equation for the evolution of baryon perturbations is given by 



^ + 2 ^-4^wW = , (30) 
ot z a at 

where in the third term in the l.h.s., the contribution from baryons has 
been dropped as it is negligible compared to the one of dark matter. 

It is convenient to define for each component the linear growth function 

D(y), 



S = D(y)S , (31) 

where y = log(a) and Sq is the initial density contrast (assuming a Gaussian 
distribution). It is also interesting to consider the so-called growth exponent 
m(y) = D (y)/D(y), where the prime denotes derivative with respect to the 
scale factor. 

Asymptotically, given that dark matter drives the evolution of the 
baryon perturbations, then they grow with the same exponent m(y). How- 
ever, their amplitudes may differ and their ratio corresponds to the so-called 
bias parameter, b = 6^/8 dm- 

It is of course, phenomenologically interesting to study the behaviour 
of 5dm, tn{y) and b as function of the scale factor a. While solving the 
differential equations for the linear perturbation, the initial conditions are 
chosen so that at a = 10~ 3 , the standard linear solution D ~ a is reached. 
In Figure 2, it is shown the linear density perturbation for dark matter, Sd m , 
as a function of a. One sees that, whereas for a = (the ACDM case), 
the perturbation stops growing at late times, for models with a > the 
perturbation starts departing from the linear behaviour around z ~ 0.25, 
the very epoch when the A term starts dominating (cf. Figure 1). In view 
of this behaviour, it is tempting to conjecture that, in our unified model, 
the interaction between dark matter and A-like dark energy is related with 
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Figure 4. Bias parameter d as a function of the scale factor, a. The solid, dotted, 
dashed and dash-dot lines correspond to a = 0,0.2,0.4,0.6, respectively. It is assumed 
that n M = 0.05, n dm0 = 0.25 and Q A o = 0.7. 



structure formation, so that for a sufficiently high density contrast (5dm » 
1), a significant energy transfer from dark matter to dark energy takes place. 
In any case, our proposal for GCG indicates that there is a connection 
between structure formation scenario and the dominance of dark energy, a 
link that ultimately results in the acceleration of the Universe expansion. 
This feature hints a possible way to understand why ttdm — just at 
recent past, the so-called Cosmic Coincidence problem. 

The behaviour of m(y) is also interesting. One can infer from Figure 3 
that from z ~ 5 to the present, the growth factor is quite sensitive to the 
value of a. For a = 0.2, m(y) increases up to 40% at present in relation to 
the ACDM case. 

Notice that m(y) governs the growth of the velocity fluctuations in the 
linear perturbation theory as the velocity divergence evolves as —HamSdm] 
it follows then that large deviations of the growth factor with changing 
a are detectable via precision measurements of large scale structure and 
associated measurements of the redshift-space power spectrum anisotropy. 

In what concerns the bias parameter, its behaviour is shown in Figure 
4. From there one can see that it also changes sharply in the recent past as 
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Figure 5. Contours for parameters b and m in the Q m -ct plane. Solid lines refer to b 
whereas dashed lines refer to m. For b, contour values are 0.98, 0.96, 0.9 from left to 
right. For m, contour values are 0.6, 0.65, 0.8 from left to right. 



a increases. This bias extends to all scales consistent with the Newtonian 
limit, hence being distinguishable from the hydrodynamical or nonlinear 
bias which takes place only for collapsed objects. Therefore, from the ob- 
servation of large scale clustering one can distinguish the non-vanishing a 
case from the a = (ACDM) case. 

The growth factor and the bias parameter at z ~ 0.15 have been recently 
determined using the 2DF survey 20 ' 21 . It is found for the rcdshift space 
distortion parameter, = 0.49 ± 0.09, and for the linear bias, b = 1.04 ± 
0.14. Notice that, as = m/b, one can obtain m = 0.51 ± 0.11. In 
Figure 5, it is shown contours for b and m in the Cl m -a plane. From the 
mentioned observational constraints on b and m, one can constrain a to 
a small but non-zero value (a ~ 0.1). However, it is important to point 
out that our study refers to the properties of the baryons whereas the 
observations concern the fraction of baryons that collapsed to form bright 
galaxies; the relation between the two is still poorly known. As far as 
parameter is concerned, one should bear in mind that this constraint is 
obtained in the context of the standard ACDM model in order to convert 
redshift to distance. Thus, a full analysis in the context of the GCG model 
is still to be performed. 

Furthermore, as can be seen from Figure 2, there is no suppression 
of 8dm at late times for any positive value of a, and hence one should not 
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expect the corresponding suppression in the power spectrum normalization, 
erg, for the total matter distribution. This was one major problem in the 
previous GCG model approach which, as pointed out in Ref. [14], cannot 
be solved even after the inclusion of baryons. In the approach developed in 
Ref. [15] and described here, one can overcome this difficulty. 

Another interesting cosmological test for our model is the study of the 
M/L ratio for clusters of galaxies. The most recent average value 22 , TO = 
0.17 ± 0.05, has been extracted from the observation of 21 clusters with 
z ~ 1. The fact that nearby cluster data seem to prefer smaller values for 
f2 m than the one obtained from WMAP data, can be interpreted as a signal 
in favour of a decaying dark matter model such as the GCG. 

4. Conclusions 

In this contribution, we have presented a setup where the GCG has been 
decomposed in two interacting components. The first one behaves as dark 
matter since it is pressure-less. The second one has an equation of state, 
px = uoxPx- It has been shown that ujx < — 1- Thus, once phantom-like 
behaviour is excluded the decomposition is unique. Apparently the model 
does not look different from the interacting quintessence models where one 
has two different interacting fluids; however, an interesting feature of our 
proposal is that it can be described through a single fluid equation. Hence, 
as far the background cosmology is concerned, we have an unified GCG 
fluid behaving as dark matter in the past and as a dark energy in the 
present. Nevertheless, when studying structure formation in this model 
one should consider it as an interacting mixture of two fluids to achieve a 
proper description. In any unified model, one expects an entangled mixture 
of interacting dark matter and dark energy. In the case of the GCG, we 
can uniquely identify the components of this mixture and the interaction. 
Moreover, we find that one does not need anything besides an evolving 
cosmological term to describe dark energy. This is consistent with recent 
studies that show that a combination of WMAP data and observations of 
high rcdshift supernovae can be described via a cosmological constant-like 
dark energy 23 . One can also consider the GCG as a decaying dark matter 
model where the decay product is a cosmological constant. 

Obviously it remains to be seen how one can obtain such a decaying 
dark matter model from a fundamental theory. Given the fact that the 
GCG equation of state arises from a generalized Born-Infeld action, it is 
possible that D-brane physics can shed some light into this issue (see eg. 
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Ref. [24]). 

Furthermore, we have demonstrated that in the context of our setup, 
the so-called dark energy dominance is related with the time when mat- 
ter fluctuations become large (6d m > 1), a possibility has actually been 
previously conjectured 25 . Moreover, we have shown that in what concerns 
structure formation, the linear regime {5dm ~ o) is valid till fairly close to 
the present, meaning that at the time structure formation begins, z c ~ 10, 
the influence of the dark energy component was negligible and that clus- 
tering occurs very much like in the CDM model. We have shown that the 
growth factor as well as the bias parameter have a noticeable dependence on 
the a parameter. We have implemented a model which exhibits a violation 
of the Equivalence Principle, as dark energy and baryons are not directly 
coupled. This may turn out to be an important observational signature of 
our approach. 
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